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Abstract 

Precision polarimetry is essential for future e'’'e“ colliders and requires Compton 
polarimeters designed for negligible statistical uncertainties. In this paper, we dis¬ 
cuss the design and construction of a quartz Cherenkov detector for such Compton 
polarimeters. The detector concept has been developed with regard to the main sys¬ 
tematic uncertainties of the polarisation measurements, namely the linearity of the 
detector response and detector alignment. Simulation studies presented here imply 
that the light yield reachable by using quartz as Cherenkov medium allows to resolve 
in the Cherenkov photon spectra individual peaks corresponding to different num¬ 
bers of Compton electrons. The benehts of the application of a detector with such 
single-peak resolution to the polarisation measurement are shown for the example 
of the upstream polarimeters foreseen at the International Linear Collider. Results 
of a hrst testbeam campaign with a four-channel prototype conhrming simulation 
predictions for single electrons are presented. 



1 Introduction 


Polarised beams are a key ingredient of the physics program of future electron-positron 
colliders (^. The precise knowledge of the beam polarisation is as important as the 
knowledge of the luminosity, since for electroweak processes, the absolute normalisation 
of expected event rates depends on both quantities to the same order. The values relevant 
for the experiments, namely the luminosity-weighted average polarisations at the interac¬ 
tion point {V)iumi, have to be determined by combining fast measurements of Compton 
polarimeters with long-term scale calibration obtained from reference processes in collision 
data. 


In particular for the International Linear Collider [^, where both beams are foreseen to be 
longitudinally polarised, it is required to control (P)iumi with permille-level precision. Two 
Compton polarimeters per beam aim to measure the longitudinal polarisation before 
and after the collision with a precision of 5VjV < 0.259fQ It should be noted, though, 
that this goal is not driven by physics requirements, but by what used to be considered 
feasible experimentally. Thus, further improvements in polarimetry would still have direct 
benefits for the physics potential of the machine. Spin tracking simulations are required 
to evaluate the effects on the mean polarisation vector caused by the beamline magnets 
between the polarimeters and the interaction point, by the detector magnets and by 
the beam-beam interaction. These effects have recently been studied in |^, concluding 
that a cross-calibration of the polarimeters to 0.1% is feasible, as well as individual 
extrapolations of upstream and downstream measurements to the e+e“ interaction point. 


The long-term average of the polarisation at the interaction point can be determined from 
reference reactions, where in particular W pair production j5]-[^ and single W produc¬ 
tion have been studied. Precisions of about 0.15% can be achieved after several years 
of data taking. These results are quite robust with regard to systematic uncertainties: 
e.g. in the measurement using W pair production, even for a conservative assumption of 
0.5% uncertainty on the selection efficiency for the signal and 5% for the background, 
the impact on the uncertainty of the polarisation measurement was found to be below 
the statistical uncertainty for an integrated luminosity of 500 fb~^ [^. Any imperfection 
in the beam helicity reversal, i.e. differences in the magnitude of the polarisation between 
measurements with left- and right-handed polarised beams, has to be corrected for based 
on the polarimeter measurements. 


The two polarimeters are located about 1.8 km upstream and 160 m downstream of the 
interaction point. Both provide non-destructive measurements of the longitudinal beam 
polarisation based on the polarisation dependence of Compton scattering and have been 
designed for operation at beam energies between 45 GeV and 500 GeV. In the order of 
10^ electron^ per bunch undergo Compton scattering with circularly polarised laser light 
which is shot under a small angle onto individual bunches. The energy spectrum of these 
scattered particles depends on the product of laser and beam polarisations as shown in 
figure and the differential rate asymmetry with respect to the sign of the circular 
laser polarisation A is directly proportional to the beam polarisation. The asymmetry 


^for typical ILC beam polarisation values of Vf.- > 80% and Vf.+ > 30% or even > 60%. 
^or positrons in case of the positron beam of the ILC which is equipped analogously. 
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expected for XV = 100% is called analysing power AV. The scattering angle in the 
laboratory frame is less than 10 /rrad, so that a magnetic chicane is employed to transform 
the energy spectrnm into a spatial distribntion. Snch a spatial distribntion is shown in 
hgure ^ for the example of the chicane foreseen for the ILC’s npstream polarimeters, 
which is sketched in hgnre 



Figure 1: Sketch of the magnetic spectrometer for the ILC npstream polarimeter 
(from [^). The Compton-scattered electrons with the lowest energies are 
deflected most by the spectrometer. 


Several options for detecting the Compton-scattered electrons behind the last magnet 
of the chicane are being considered. The baseline solntion is a gas Cherenkov detector 
consisting of 20 channels, each 1 cm wide. The rate asymmetry Ai with respect to the 
laser helicity is determined in each channel i as 


Ai 


ATf+ iV^’ 


( 1 ) 


where and AT" are the connt rates in channel i for right- and left-handed laser po¬ 
larisation, respectively. The connt rate and asymmetry for each channel for the baseline 
solntion as well as a detector design with a three times larger channel number are shown 
in hgure and hgure respectively. 


Per-channel measurements of the beam polarisation are then obtained by comparison with 
the predicted analysing power AVi of the corresponding channel. The predictions of AVt 
are based on the well calculable cross-section for Compton scattering]^ plus a Monte-Carlo 
modelling of the magnetic chicane and the detector geometry and response. The resulting 
per-channel polarisation measurements are then combined in a weighted average, where 
the largest weights Wi are given to the channels with the smallest statistical uncertainties. 


which is equivalent to the channels with the largest analysing power 11 


Wi = 


{V ■ AV, 


1-2 


^With radiative corrections of less than 0.1% currently calculated to the order of 


10 


( 2 ) 
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(a) 



(b) 



Figure 2: (j^ Differential cross-section in dependence of the energy of the scattered elec¬ 
tron for XV = ±80 %. (|^ Spatial distribution of the Compton scattered elec¬ 
trons at the detector location behind the magnetic chicane for XV = ±80%. 
The hlled area indicates the resulting asymmetry (secondary y-axis). The bot¬ 
tom plots show the resulting mean number of Compton electrons per detector 
channel for XV = ±80% and the corresponding asymmetry assuming (|^ 20 
channels with a width of 10 mm and ([^ 60 channels with a width of 3 mm. 


This measurement principle has already been used successfully at the SLC polarimeter, 
which was the up to now most precise Compton polarimeter, reaching 5V jV = 0.5% 

The gas Cherenkov approach is robust against high data rates, a harsh radiation en¬ 
vironment, and possible backgrounds. The basic detection mechanism is intrinsically 
proportional to the number of electrons passing a single channel simultaneously, since 
the amount of emitted Cherenkov radiation does not depend on the energy of the par¬ 
ticles once they are relativistic. The main source of non-linear behaviour are thus the 
photodetectors. Table [T] shows the anticipated uncertainty budget for the ILC Comp¬ 
ton polarimeters assuming a Cherenkov detector, and compares it to the uncertainties 
achieved at the SLC polarimeter. 
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The most challenging improvement has to be achieved with respect to the detector lin¬ 
earity. In order to reach the goal of not more than 0.1% contribution from this source to 
5VIV, non-linearities have to be monitored and corrected for at the level of a few per- 
mille [^. A corresponding calibration system has been developed recently [^, matching 
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source of uncertainty 5V jV 



SLC achieved 

ILC goals 

laser polarisation 

0.1% 

0.1% 

detector alignment 

0.4% 

0.15%-0.2% 

detector linearity 

0.2% 

0.1% 

electronic noise and beam jitter 

0.2% 

0.05% 

Total 

0.5% 

0.25% 


Table 1: Uncertainty goals for the polarisation measurement at the ILC [4|. For com¬ 
parison, the systematic uncertainties determined for the SLC polarimeter are 
also given 
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a prototype of a gas Cherenkov detector which has been successfully operated in testbeam. 
With this prototype, the alignment requirement has nearly been reached and electronic 
noise was found to be at a negligible level [^. The influence of beam jitter will also be 
much smaller in the ILC case, since the whole accelerator has been designed to limit 
the transverse jitter to smaller than 10% of the beam sizes in order to achieve the 
design luminosities. Thus, it seems feasible to reach the envisaged precision goal with the 
baseline design. 


The effect of non-linearities would be eliminated completely if the Compton-scattered 
electrons could be detected individually, e.g. by a Silicon pixel detector, which would also 
increase signihcantly the possibilities to control the alignment. Here however, R&D is still 
required to reach a proof-of-principle level, where the main concern is the high local data 
rate, followed by radiation issues. Employing pixel detectors would enable in addition the 


measurement of transverse beam polarisation 15 


An alternative Cherenkov detector concept which promises a signihcant reduction of the 
effects of non-linearities and misalignments is the subject of this paper. It is organised as 
follows: In section]^ the concept of polarimetry with Compton electron counting is intro¬ 
duced and its benehts are demonstrated using the example of the upstream polarimeter 
of the International Linear Collider. Section discusses the design of a suitable detector 
matching the requirements and introduces a corresponding prototype. In section we 
present results from operating this prototype in testbeam before concluding in section 


2 Polarimetry with Compton electron counting 


The detector concept discussed in this paper aims to resolve individual Compton electron 
peaks in the detected and digitised spectrum of the Cherenkov light. This would allow to 
build a “self-calibrating “ detector: by determining the distance between these peaks, the 
response to single Compton electrons including the amount of detected Cherenkov light as 
well as the gain of the photodetector could be monitored in-situ. Even better, the shape of 
such a spectrum could be exploited to directly determine the number of Compton electrons 
in this channel for the calculation of the polarisation asymmetry, without the need to 
calibrate the photomultiplier gain, its quantum efficiency and geometrical acceptance for 
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the produced Cherenkov light. It should be noted that this goes beyond the single-peak 
resolution offered for low count-rates by Silicon photomultipliers (SiPMs), where each peak 
corresponds to a fixed number of detected photons. Thus, the single peaks of a SiPM do 
not give direct information on the Compton spectrum, but would need to be unfolded 
for the statistical and wavelength-dependent effects of the Cherenkov light production, 
propagation and detection. When operated at the higher count-rates required to resolve 
individual Compton electron peaks, SiPMs exhibit a significantly higher intrinsic non¬ 
linearity and temperature dependence compared to conventional photomultiplier tubes. 
Therefore, we do not consider SiPMs further in this concept study. 


input parameters: 


simulated quantity: 


method: 


mean number of 
Compton electrons 

mean number of photo¬ 
electrons per Compton e“ 

gain, gain fluctuations 
non-linearity parameters 


QDC resolution 



Poisson distributed 

Poisson distributed 

Gaussian distributed 
NL transfer function 

histogram binning 


Figure 3: Schematic overview of the Monte Carlo simulation predicting the digitised 
detector response, with the input parameters for each simulation step in the 
left column and the method to implement the corresponding statistical effects 
in the right column. 


A dedicated Monte-Carlo simulation has been developed to predict the digitised photomul¬ 
tiplier response for a given average number of Compton electrons per detector channel. As 
illustrated in figure]^ the simulation starts by determining the actual number of Compton 
electrons per channel based on a Poissonian distribution around the mean value. For each 
Compton electron, the actual number of photoelectrons created and emitted from the 
photocathode of a photomultiplier is determined, again assuming Poissonian statistic^ 
The charge amplification process is then modelled including Gaussian gain fluctuations 
(in the following also referred to as “noise “) and, optionally, a non-linear transfer func¬ 
tion. Finally, the resulting charge is digitised with a fixed resolution. Any potential 
imperfections of the digitisation step are implicitly included in the non-linearities of the 
amplification step, since they have been found to be tiny for our actual digitiser 16 


Further noise from the read-out electronics is assumed to be negligible in the scope of the 
detector design studies presented here, based on the testbeam experience with the gas 
Cherenkov prototype 1^. 


^This is equivalent to a Poissonian (or for large numbers Gaussian) distribution of produced Cherenkov 
photons and binomial statistics for the detection probability. 
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2.1 The case for single-peak resolution 

In this section, we will first derive the conditions reqnired to resolve individnal peaks in 
the digitised photomnltiplier spectrnm. We will then show how the single-peak resolntion 
can be exploited for a reliable estimate of the mean nnmber of Compton electrons per 
channel, even if the detector behaviour is not perfectly linear. 


2.1.1 Requirements for single-peak resolution 

Two factors are important to resolve individual Compton electron peaks: the average 
number of Compton electrons Nc.e. per detector channel, and the number of photoelec¬ 
trons per electron p := Np^eX^c.e. = !)• To illustrate this, £gure|^ shows a simulation of 
digitised spectra for Nc.e. = 5, with a yield of p = 7 and p = 300, respectively. Here, a 
gain g of 4 ■ 10^ with a hxed noise level of ^ = 1 % followed by charge digitisatior0 with 
a resolution of wqdc =200 fC is assumed as example. For the conhguration with p = 7, 
the Poissonian distribution of actual electrons per bunch crossing is smeared into one 
broad peak. For p = 300, however, peaks for the individual actual numbers of Compton 
electrons ic.e. can be easily distinguished. 




Figure 4: Simulated spectra for 5 Compton electrons in the detector with (j^ 7 and (j^ 
300 detected photons per Compton electron. The simulation assumes that the 
photodetector signal is amplihed with a gain of 4 ■ 10^ with a hxed noise level 
of 1 %, followed by digitisation with a QDC with a resolution of 200 fC. 

To resolve individual peaks, the separation s between them should be larger than the 
peak width by a factor k, which can be phrased as a requirement of 

s > FWHM = fc ■ a, (3) 

with a the width of a Gaussian approximation of an individual peak. 

®In the following we will use the abbreviation QDC (charge sensitive analogue-to-digital converter). 
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The separation s is determined by p, the photomultiplier gain g and the digitiser resolution 


wqbc as illustrated in hgure 5a 


s = p ■ 


9 • e 

■W^QDC 


= p-Q, 


( 4 ) 


where e is the electron charge. The peak width is given by a convolution of the fluctuations 
in the actual number of photoelectrons per Compton electron around p, the noise Ag/g 
of the photomultiplier, and Q as dehned by equation as 


a = 


ic.e. ■ P + Q^ic.e. ■ P 




( 5 ) 


The hrst two terms are attributed to the statistics of the photoelectron creation and their 
amplihcation, respectively. The term 1/12 accounts for the digitisation error, which is 
negligible for any of the conhgurations considered here. To meet the separation require¬ 
ment given by equation for at least all peaks up to the mean number of Compton 
electrons Nc.e., a combination with equations [5|4| shows that p needs to fulhl 


I'Ve... + 

P > - - -^- A - 

i - {Nc... ■ fr 

With a choice of fc = 2.35 for a separation by at least the full width at half maximum, 
at least p > 30 is required to resolve the case Nc.e. = 5 for the previous example from 
hgure 1^ For the case of p = 300 as assumed in hgure [4b| up to Nc.e. = 29 can be resolved. 

Obviously, the modelling of the photomultiplier noise is important for a correct estimation 
of the peak separation capability for a given setup. While the example in hgure [^assumed 
a constant noise level of 1%, a more realistic approach based on a statistical modelling of 
the amplihcation process jrn leads to a dependency of ^ on the number of photoelectrons 
which are produced at the cathode. This is taken into account in the simulation studies 
presented below and results in changes of about 10% in the value of p required to resolve a 
certain Nc.e.- Thus equationwith a constant noise level remains a useful approximation. 


( 6 ) 


2.1.2 Reconstruction of A^"c.e. 


The Monte-Carlo simulation described above has been employed to obtain QDC spectra 
for various conditions. As example, the histogram in hgure ^ shows the predicted QDC 
spectrum for A^c.e. = 20 assuming p = 300, g = 2.4 ■ 10® with detailed modelling of the 
noise, perfect linearity and wqdc = 200 fC. There are two possibilities to reconstruct the 
average number of Compton electrons from this spectrum: either by taking the mean (or 
maximum) of the whole distribution as in the case with no single-peak resolution, or by 
extracting A"c.e. from a ht to the detailed shape of the spectrum. 


A suitable ht function is the sum of Gaussians with mean q 2 -ip^^^ and width g2-ipeak+i) 
where each Gaussian describes one of Speaks peaks in the spectrum. Since the number of 
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Figure 5: (j^ Illustration of peak position, width and separation, for the hrst three peaks 
in a simulated QDC spectrum. ^ QDC spectrum for Nc.e. = 20 and p = 300 
(black histogram) htted according to equation(blue curve). 


electrons in a channel is expected to follow a Poissonian distribution, the height of the 
peaks is set to the expectation for ipeak electrons from a Poissonian distribution of Nc.e. 
electrons and scaled to the number of events Nevents in the spectrum: 

“Speaks 

f{xQDc) = Nevents ' ^ Pois(ipeak, Nc.e) ' Gaus(a:Qi5c, 92.ipeak) ?2.ipeak+l)- (7) 

^peak~0 


The free parameters of the £t are the mean and width of all Upeaks peaks and the central 
value Nc.e. of the Poissonian. While studying the distance between the mean values 
of the individual peaks gives access to the gain linearity, the central value Nc.e. of the 
Poissonian is the parameter of interest for the polarisation measurement. An example for 
a £t of equation to a spectrum simulated with Nc.e. = 20 is shown in hgure ^ for 10^ 
accumulated individual measurements. 


The number of parameters to be htted grows with the number of electrons per detector 
channel. For Nc.e. = 20, peaks up to ic.e.~40 need to be described in the QDC spectrum 
and consequently a £t function with ?«80 parameters is required. For such a large number 
of free parameters, a careful choice of start values for the £t is essential. Since the distance 
between neighbouring peaks is constant for a perfectly linear photomultiplier and changes 
only gradually for non-linear gain, an initial estimate for s can be obtained from the 
discrete Fourier transform of the spectrum. From the separation, start values for the 
Gaussian parameters can be determined by calculating the expected mean and width. 
Since s corresponds to the scale factor between initial Gompton electrons and QDG bins, 
a start value for Nc.e. can be obtained by fitting a Poissonian P{x') with x' = 
the spectrum. 


Figure 6a shows the relative deviation of the fit result for Nc.e. from its true value for 
up to Nc.e. = 35 in case of a perfectly linear response (blue) and in the presence of a 
photodetector non-linearity of 0.5 %. 
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Figure 6: Relative deviation of the measured Nc.e. from its true value when Q exploit¬ 
ing single peak resolution and (|^ when relying on the mean of the QDC 
spectrum, respectively. The central lines are the mean deviations for 50 simu¬ 
lation runs per value of Nc.e. with p = 300, the hlled bands represent the RMS 
around the central value. In case of a perfectly linear detector response, both 
methods perform equally well (blue). When a photodetector non-linearity of 
0.5% is simulated (using non-linearity functions with different randomly se¬ 
lected parameters for each run), the single-peak based reconstruction shows a 
superior performance (orange). 


Up to ~25 Compton electrons, the added non-linearity has little effect. For larger num¬ 
ber of electrons, the description of the spectrum by the £t start values does not model 
the spectrum shape sufficiently well anymore and the fit results starts to degrade. For 
comparison, figure 6b shows how much a photodetector non-linearity of 0.5 % affects the 
conventional method using the mean of the QDC spectrum, which is the only method for 
the calculation of the polarisation asymmetry when no single-peak resolution is possible. 


2.2 Application to polarisation measurements at the ILC 


The results in £gure|^ where based on accumulating 10^ individual measurements. With 
the baseline beam parameters of 1312 bunches per bunch train at 5 Hz bunch crossing 
rate, 3.936 ■ 10® measurements could be collected in 10 minutes, providing a large enough 
data sample so that the fit performance is not limited by statistics. Since the non-linearity 
is not expected to change rapidly, the actual polarisation measurement can be done on 
much smaller datasets, with a determination of the peak positions from the data taken in 
the previous 10 minutes or, in an offline analysis, with a moving average over data taken 
before and after the individual short-time dataset. 


To evaluate the benehts of the £t procedure described above for polarisation measurements 
at the ILC, a detector array of 60 channels with a width of 3 mm and 0.33 mm inter¬ 
channel spacing was simulated. The fast Linear Collider Polarimeter Simulation LCPolMC 
is used [4,13 to generate Compton events and subsequently track the Compton electrons 
through the polarimeter chicane to the detector. The beam and laser parameters used 
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Figure 7: Polarisation calculated for each detector channel i for a simulation of 80% po¬ 
larisation with Q 0.2% and (|^ 1.6% photodetector non-linearity, respectively. 
The orange markers show the polarisation Vc.e.,i calculated from the primary 
Compton electrons, the blue markers the calculation 'Pqdc,! from the mean of 
the QDC spectrum, and the green markers the results Vm,! for the number of 
Compton electrons determined from the £t to the QDC spectra. The fluctu¬ 
ations around channel 30 are due to the zero-crossing of the asymmetry (c.f. 
hgure 2d). 


for the simulation are chosen according to the expectations at the polarimeter locations 
for the ILC TDR beam parameters [^. The resulting number of Compton electrons per 
channel is displayed in figure It is below Nc.e. = 30 for all channels, and consequently 
the requirement for single-peak resolution given by equation is fulhlled. 


2.2.1 Robustness against non-linear detector response 


Based on figure [Mf the same simulation chain as before is applied to all channels for both 
laser helicities and the average number of Compton electrons per channel is extracted 
both from single-peak htting and from the mean of the QDC spectrum. The polarisation 
is then reconstructed as described in section [U 


Figure 7a compares the polarisation measurements obtained via the two reconstruction 
methods in case of a rather small non-linearity of 0.2% with the ideal result obtained 
directly from the true Compton electron spectrum. While the results from the three 
different methods exhibit no striking difference in this case, a larger non-linearity of 1.6 % 
as shown in hgure 7h distorts the QDC mean and therefore the measured polarisation 
significantly, with a relative deviation of 0.39 %, while the polarisation calculated from 
the fit results appears to be unaffected, with a deviation of 0.01 %. 


For a more systematic survey of the impact of non-linearities, the same procedure as 
for figure Ta was repeated for non-linearities from 0% to 4% in 0.2% steps. For each 
step, random parameters for 100 different functions for the photomultiplier non-linearity 
were picked and used for the spectrum generation. The mean and RMS of the deviation 
between the polarisation determined from the mean of the generated QDC with respect to 


11 













Figure 8: Deviation between the resnlt from the polarisation Vc.e. calcnlated from the 
primary Compton electrons and 0 the polarisation Vqbc determined from 
the mean of the QDC spectrum, and (|^ the polarisation "Pfit calculated with 
the number of Compton electrons determined from the £t of the QDC spectra, 
using all detector channels (blue, solid line) or only those channels i whose 
polarisations Pfit,i agree within ±1%-Pfit with the polarisation calculated from 
all channels (green, dash-dotted line). 


the result for the primary electron calculation is shown in figure The error budget for 
the detector linearity contribution to the polarisation measurement is 0.1 % (see table [^. 
The method using the mean of the QDC surpasses the allocated limit for non-linearities 
> 0.4%. The results for the polarisation calculation from the htted number of Compton 


electrons is shown in hgure 8b Fitting the QDC spectra instead of using the mean can 


compensate the photodetector non-linearity well enough to stay within the error budget 
for non-linearities up to ~ 2.2 %. 


This can be further improved by introducing a second step in the calculation of the 
polarisation to eliminate the contribution from outlier channels with poor £t performance. 
A rough hrst approach to demonstrate this is a recalculation of the polarisation using 
only the channels which agree within ±1 % • Pat with the polarisation calculated from all 
channels. The fraction of channels consequently used in the calculation is larger than 90 %. 
Using less channels will decrease the statistical precision slightly. However, the benehts 
of this approach in dealing with larger detector non-linearities, compared to the use of 
all channels, are clearly visible in hgure where the deviation between the polarisation 
thus obtained with respect to the result for the primary Compton electrons is shown in 
green. Even for photodetector non-linearities of 4%, this deviation is below 0.03%, i.e. 
well within the error budget of 0.1 %, which would be a contribution towards reducing the 
overall systematic uncertainty on the polarisation measurement. For the method using the 
QDC mean, such an approach would not lead to such improvements, since all channels are 
systematically affected by the non-linearity, as one can see when looking at the example 
in hgure 
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channel 


(b) 


Figure 9: (j^ Deviation of the polarisation obtainable from the primary Compton elec¬ 
trons in dependence of detector rotation in the horizontal plane for a simula¬ 
tion of 80 % polarisation, and 0 the polarisation calculated for each detector 
channel i without detector misalignments and in case of a 1.4mrad tilt in the 
horizontal plane. 


2.2.2 Detector alignment 


Next to non-linearities, the detector alignment with respect to the Compton electron 
fan is the other large contribution to systematic uncertainties on the polarisation mea¬ 
surements that should be taken into account in the design of the polarimeter detectors, 
c.f. table The largest contribution originates from the alignment in the plane of the 
magnetic chicane, since it directly affects the analysing power predicted for each channel. 


Rotations of the detector in the deflection plane of the chicane (i.e. in the horizontal 
plane in case of the upstream polarimeter) affect the path of the electrons through the 
channels. The amount of Cherenkov light per electron and its propagation to the pho¬ 
todetector could be altered, thus changing p. This effect has to be found negligible for all 
simulated angles up to 5°. More importantly, such horizontal tilts can introduce cross-talk 
if a Compton electron passes through more than one channel. Figure |9^ shows the impact 
on the polarisation: tilts of the detector in the horizontal plane result in deviations of 
0.07%/mrad. For the ILC polarimeters, a tilt alignment of 1 mrad or better is consid¬ 
ered feasible, thus the corresponding systematic uncertainty is within the allocated error 
budget. Horizontal rotations of the detector have the strongest impact on the channels 
with the smallest asymmetries near the zero-crossing of the asymmetry, as can be seen 
in hgure This characteristic pattern can be resolved with high granularity as in case 


of the 60 channels considered here. Consequently, the systematic uncertainties could be 
further reduced by identifying the presence of misalignments and correcting for them or 
excluding the most strongly affected channels. 


Misalignments of the detector position in the horizontal direction can be controlled by 
monitoring the fraction of Compton electrons in the detector channel at the Compton 
edge. Figure 10a shows how the fraction of Compton electrons in the detector channel 
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Figure 10: (j^ Fraction of Compton electrons expected in the detector channel contain¬ 
ing the Compton edge as a function of a horizontal offset for "P = 30 % and 
80 % polarisation, averaged over the measurements with both laser helicities. 
A hrst order polynomial has been htted to the 80 % case. e Polarisation 
"^offset,! calculated for each detector channel i for a simulation of 80 % polar¬ 
isation with a horizontal misalignment of ~ 50 /im (green). For comparison, 
the blue (orange) markers show the polarisation Pat,! (Pc.e.,i) obtained with¬ 
out misalignment from the £t to the QDC spectra (the primary Compton 
electrons). 




— Vi from C.e. 

Vi from fit 
_ Vi from 


P,., =(80.007±0.003)% 
Pfi. =(80.009±0.004)% 
^.f<^.= C79.979±0.004)% 


containing the Compton edge over the total number of detected electrons changes as a 
function of a horizontal offset of the detector. 

_ -^C.e(pdge) 

^ “ Auchan 


Averaging over the measurements with both laser helicities eliminates any dependence 
on the polarisation, as can be seen in hgure 10a A £t of a linear function to this ratio 
yields a slope m of —0.343 It should be noted that the expected number of Compton 
electrons hitting the front face of each detector channel for a given offset is independent 
from intrinsic detector uncertainties and thus the precision of m is only limited by the 
knowledge of the held of the chicane magnets, which is expected to be much better than 
the detector-beamline alignment. 


In case of single-peak resolution, ANc e{i) = 0.1 % is estimated for the determination of 
the Compton electron number per channel based on the performance shown in figure 
which holds even in case of a not perfectly linear detector response By propagating the 
assumed uncertainty for each channel to the ratio y, Ay = 7.9- 10“^ % is obtained, which 
translates into a horizontal alignment precision of Ax = \Ay/m\ = 2.3/rm. 


Such an alignment is sufficient to meet the requirements on the contributions to the 
systematic uncertainty of the polarisation measurement from the detector alignment. 


^The same method can also be applied when relying on the mean of QDC spectrum only. In that 


case, however, the precision degrades significantly in case of non-linearities, c.f. figure 6b 
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Figure 10b shows the impact that of horizontal misalignment of 50 /im in a simulated 
polarisation measurement. The deviation of AV/V = 0.03% caused by the offset is well 
within the allocated error budget. Thus htting the single-electron-peak spectra is also a 
promising approach towards meeting the requirements for the detector alignment. 


3 Design of a quartz Cherenkov detector 

For the baseline gas Cherenkov detector with 1 cm wide channels, the mean number 
of Compton electrons per channel reaches up to A'^c.e. — 100, while the average yield 
of photoelectrons per Compton electron is p = 6.5 [^. This conhguration is far from 
fulhlling equation In order to develop a Cherenkov detector suitable to achieve single¬ 
peak resolution, there are two possibilities: The number of Compton electrons per channel 
can be reduced by a factor 2-3 by building smaller channels for the polarimeter. In 
addition, it is mandatory to increase the number of photoelectrons per Compton electron 
substantially. As an approach to achieve the latter, the use of quartz as Cherenkov 
material is considered. The higher refractive index compared to Cherenkov gases results 
in a much higher light yield, as described by the Frank-Tamm formula 18 . Compared 


to e.g. perfluorobutane gas, the use of quartz would translate into an approximately 200 
times higher intensity of the emitted Cherenkov light for relativistic particles. However, 
this will be somewhat mitigated by other effects, such as increased absorption as the 
produced light travels to the photodetector. In this section, we discuss a simulation study 
of these effects and derive a suitable design for a prototype detector. 


3.1 Geant4 simulation 


To study the feasibility of a quartz detector for ILC polarimetry and investigate different 
design options, as well as for comparison to data from a prototype detector, a detailed 
simulation of the detector concept based on Geant4 19 has been developed. The 


simulated geometry is depicted in hgure 11 The main elements of the simulated setup 


are quartz blocks surrounded by a thin layer of air and/or aluminium foil, the entrance 
window and photocathode of a photomultiplier tube (PMT) and a layer of optical grease 
between the quartz block and the PMT window. To increase the space available for the 
photodetectors and their readout, every other quartz block is flipped with respect to its 
neighbours. 

To produce Cherenkov light, electrons are shot through the quartz. For the Cherenkov 
photons, all relevant processes are simulated, in particular absorption according to the 
absorption length of the relevant material, Rayleigh scattering, and boundary processes at 
the surface between two different media. For the detector material, the optical properties 
of Spectrosil® 2000 20 have been implemented, and for the optical grease the proper¬ 
ties of Cargille fused silica matching liquid code 06350 and 50350 21 . To describe the 


behaviour of the photons at boundaries between materials, the UNIFIED model 
Geant 4 is used. 


22 of 


Various geometrical properties were varied to study their impact on the light yield: 
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Figure 11: Sketch of the simulated geometry. The left side Q shows some components 
of the simulation, in addition demonstrating the added space for electronics 
due to the rotation of every second quartz block. For a better overview, on 
the right side a the same geometry is depicted in side- and top-view, along 
with the geometrical parameters varied in the design studies: the channel 
dimensions width Lx, height Ly, length L^, the incidence angle ax', as well as 
the space between quartz and electrons dquartz, between electrons and PMT 
dpMT, and the distance between the end of the channel and the point where 
the electron crosses the detector axis dend- 


• A larger channel width Lx reduces the number of reflections at the channel walls a 
photon undergoes before detection, which increases the probability to arrive at the 
photodetector. 

• A larger channel height Ly increases the path length of the electron crossing the 
quartz channel and consequently the amount of Cherenkov light produced. Addi¬ 
tionally, as in the case of larger width, the number of reflections is reduced. Both 
effects result in an enhanced photon yield. 

• A shorter channel length L^ reduces the distance the photons have to travel before 
detection, which increases the light yield by causing less photons to be absorbed 
inside the quartz. 

• The incidence angle ax of the electron into the quartz channel also changes the 
path length of the electron crossing the quartz channel for light production, with a 
higher light yield for small angles. Another, though less pronounced, impact on the 
number of photons is due to the fact that for angles close to the Cherenkov angle, 
a larger fraction of the photons can reach the detector without being reflected on 
the narrow side faces. The spatial distribution of the light on the photodetector 
surface is the most uniform under such angles, whereas for smaller or larger angles 
the intensity distribution is less homogeneous, as illustrated in figure [T^ This effect 
could possibly aid in the angular alignment of the detector if a photodetector with 
multiple anodes per quartz channel was used. 
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Figure 12: Light distribution on the photocathode for different incidence angles of the 
electron. The position on the photocathode surface is coloured according to 
the relative light concentration: the regions with most light for each angle 
are shown in yellow, with a colour gradient to the regions with the least light 
in dark blue. 


When combining the parameters listed above to chose a suitable detector geometry, some 
considerations have to be taken into account in addition to the light yield. Small channel 
width lead to less Compton electrons per detector channel and therefore aid in fulhlling the 
condition for reaching single peak resolution. Another consideration is that long, narrow 
channels and large angles move the photodetector further from the Compton electrons 
to protect them from contact with high energetic particles. To take this into account, 
the channel height and the electron entrance point were varied along with the length, 
such that the space between the lower end of the quartz and the electron plane stayed 
constant at dquartz = 10.0 mm and the shortest distance between the electrons and the 
photomultiplier was dpMT = 30.0 mm. This translates into a channel height Ly of 

7- _ 2 ■ L2: ■ sin^Oa;) dquartz 

ly / \ 

2 ■ COS(Q:a;) 


( 9 ) 


and a distance dend between the end of the detector channel and the electron crossing of 
the detector axis of 


dpnd — 


dquartz “1“ L/y * COs(o^ 3 ^) 


( 10 ) 


for each simulated channel length and incidence angle a^- The results for channels 
with a width of = 5 mm considering different lengths in the range 50 mm - 300 mm and 
angles in the range 30°- 60° are displayed in figure In this range of parameters, the 
best photon yield of ~3000 photons was achieved for ax = 60° and Lz = 160.0 mm, which 
allows a channel height of Ly = 197.2 mm while keeping the chosen safety margin between 
incoming electrons and readout electronics. A mechanical setup with such nearly square¬ 
shaped channels rather than elongated ones might pose some challenges, e.g. covering the 
readout face of high but narrow with conventional PMTs does not seem reasonable with 
standard photomultiplier geometries. However figure [T^ shows that there is a large variety 
of configurations with less extreme aspect ratios which still offer an impressive light yield 
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Figure 13: Photon hits for different quartz lengths, heights and angles relative to the 
electron chosen such that conditions |9] and [To] are fulhlled. The dashed blue 
lines indicate the length and height used for the prototype detector. 


in the order of ~1500 to ~2000 photons and which would be more straightforward to 
implement. 

For the construction of the prototype detector (section [3^ , smaller dimensions of Ly = 
18.0 mm and = 100.0 mm were chosen, despite the prediction of less (~600) photon 
hits. This choice was partly made to maintain the angular flexibility to scan the full range 
from 30°- 60° to conhrm the incidence angle’s effect on the light distributions, but also to 
match the dimensions of the available photodetectors, since the time and cost for tailor- 
made photomultipliers were not considered reasonable before a successful demonstration 
of the concept. 


The photon hits given above are the number of photons that reach the photodetector 
surface. Most common photodetectors have a detection efficiency well below 100 %. A 
more realistic estimate of the photon yield which could be achieved will therefore have 
to take the characteristics of the photodetector into account to evaluate the number of 
detected photons. Figure 14a shows the quantum efficiencies of two PMTs with good 
sensitivity in the UV range. 


The amount of measurable light can be further influenced by the choice of the other ma¬ 
terials used in the detector, such as the type of optical grease or the surface finish of 
the quartz bars. The impact that the choice of photodetector and materials has on the 
amount of measurable light is illustrated by figure 14b Despite the fact that layer of 


grease is only 1 mm thick, the UV component of the Cherenkov light is severely affected 
by its short absorption length. This clearly indicates that the choice of a suitable opti¬ 
cal grease to couple the photodetectors to the quartz is also of great importance. The 
application of another optical grease with lower absorption in the UV range would allow 
to increase the amount of measurable light for a chosen geometry. To a smaller extent, 
further improvements could be achieved by polishing the quartz surface. Table lists 
the predicted light yield for both the geometry with maximum predicted light yield from 
figure 13 as well as for the prototype discussed section 3.2 (with a size of 5 mm x 18 mm 
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Figure 14: Effect of photomultiplier choice: the left plot (|^ shows the quantum effi¬ 
ciencies of two UV sensitive PMTs [^. On the right (|^, the wavelength 
spectrum of the Cherenkov photons and the effects of absorption inside the 
detector and the different quantum efficiencies is shown. 


simulated geometry: 


optimised prototype 

ground polished ground polished 


photon hits 

3344 ± 58 

3876 ± 62 

576 ± 24 

651 ± 26 

detected photons 

322 ± 18 

373 ± 19 

57 ±7 

60 ±8 


Table 2: Light yields (with statistical errors) for an optimised geometry and the pro¬ 
totype detector geometry. The values for “photon hits“ refer to the number 
of photons reaching the readout at the end of the channel per incident elec¬ 
tron, while “detected photons “ are those remaining after absorption in Cargille 
50350 optical grease and application of the quantum efficiency of a Hamamatsu 
R7600U-03 PMT. 


X 100mm and for an incidence angle of ax = 45°). Despite the fact that only 0(10%) 
of the photons are detected with the material choices considered in these hrst simulation 
studies, the predicted yield of 60 - 300 detected photons per incident electron are promis¬ 
ing for achieving a single Compton electron resolution as discussed in section For the 
prototype geometry along with a PMT gain of g = 4- 10^ and a digitiser with a resolution 
of 200 fC, this would correspond to an allowed range of up to Nc.e. = 10 electrons per 
channel, while a light yield of 200 - 300 as predicted for other geometries would allow 
Nc.e. =28 - 42 electrons per channel. 


3.2 Construction of the prototype 


The simulation studies indicate that the desired photon yield can be reached. To verify 
the simulation’s predictions, a prototype detector with adjustable incidence angle was 
built and operated at the DESY 11 testbeam. 


The prototype detector consists of four quartz bars, with every second quartz block flipped 
with respect to its neighbour, as described in section |3.1[ The synthetic fused silica 
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brand Spectrosil® 2000 20 was selected as the most suitable material, based on the high 


radiation tolerance compared to other natural or synthetic fused silica brands 24 , the low 


number of optical impurities and the high transmission in the UV range. The dimensions 
of each quartz block were chosen to be 5 mm x 18 mm x 100 mm, which is a suitable 
compromise between a high light yield and off-the-shelf photodetector geometries. Each 
quartz block was wrapped in aluminium foil as shown in hgure |15a 


The photodetectors for the prototype detector were chosen according to two main criteria: 
high quantum efficiency at short wavelength, where the Cherenkov light intensity is the 
highest, and the dimensions of the sensitive area. Square photomultipliers with an active 
area of 18 mm x 18 mm, with a four-anode readout (Hamamatsu R7600U-03-M4 and 
R7600U-04-M4) were employed to read out the two quartz blocks pointing into the same 
direction, with two anodes covering one quartz block as illustrated in figure 15b This 


layout allows to have the four detector channels directly side-by-side. Figure [T^ shows 
two quartz blocks and one of the PMTs compared to a 2-Euro coin. 

To limit the occurrence of total reflection at the boundary between quartz bar and PMT 
window, they need to be coupled with optical grease. For the testbeam campaign, Cargille 
06350 was used. 

The mechanical setup to hold the quartz bars and photomultipliers in place with ad¬ 
justable angle between them is shown in hgure 16a The angle between the upper and 


lower channels can be controlled by the movement of a stepping motor with a base step 
angle of 0.9°, which is transformed to 0.5° steps by a dedicated gear mechanism contained 
in the mechanical setup. 



Figure 15: (j^ Quartz blocks with the coupling piece used to attach them to the mechan¬ 
ical holder, in the process of being wrapped in aluminium foil. (|^ Sketch 
to illustrate how one photomultiplier with four readout-anodes can be posi¬ 
tioned to provide two measurements for two quartz channels each. ([^ Quartz 
blocks and a photomultiplier with a coin for size comparison. 
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(a) (b) 


Figure 16: In the technical drawing (|^, the positioning of the four quartz blocks and 
two photomultipliers used in the prototype setup are marked. On the photo 
of the testbeam setup (|^ a sketch of the beam path is overlaid on a close-up 
of the light-tight box with the detector inside on the x-y-table and two of 
the trigger scintillators. 


Single electron measurements at the DESY test- 
beam 


The quartz prototype detector was characterised with single electrons at the DESY test- 
beam 25 . The electrons are produced by inserting a metal conversion target into a 


bremsstrahlung beam from the DESY II synchrotron. As a converter target for the gen¬ 
eration of the electrons, a copper foil of 5 mm thickness was chosen. With a dipole magnet 
and a collimator electrons with an energy of 3.75 GeV were selected to achieve a high rate 
for the data taking. In combination with an area collimator with aperture of 5 x 5 mm 
located 5 m upstream of the prototype detector, this resulted in single electrons reaching 
the prototype detector with a frequency of «ilkHz - 2 kHz. 

The quartz detector was placed inside an aluminium box with foam rubber seals to shield it 
from external light and placed on a base plate movable by two linear precision translation 


stages, as shown in hgure 16b 


The main components of the setup and the readout chain during the testbeam campaign 
are sketched in hgure |17a[ The photomultipliers are operated with a supply voltage of 


900 V, leading to amplihcation with a gain of ~ 1.5 ■ 10® (c.f. section 4.5). The charge 


signals of the photomultipliers are read out with a 12-bit charge-sensitive analogue-to- 
digital converter (QDC) which offers two resolutions, 25 fC and 200 fC. The digitisation 
length and cycle are steered by an external gate. To trigger on beam electrons, the 
generation of the gate signal was started in case of a four-fold coincidence between two 
trigger scintillators before and after the detector, respectively. A total number of eight 
QDC channels was used to digitise the signal from the two four-anode photomultipliers. 
The conventions for channel naming and anode numbering are illustrated in hgure |17b 


An example for the QDC signal recorded for a PMT anode is shown in 18a 
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(b) 


Figure 17: (|^ Experimental setup at the testbeam. The main components of the setup 
and the path of the readout signal path are sketched. B Channel names 
and anode numbering used for the description of the testbeam data. 


As a hrst step the detector response and alignment of the testbeam setup were studied, 
which consequently allowed to compare the measured data to simulation and evaluate the 
light yield. 


4.1 Detector response 


The spectrum of each photomultiplier anode contains events in which the data acquisition 
was triggered, but no electron entered the respective channel. Therefore, in addition to 
the Cherenkov signal another peak is present at the low edge of the spectrum, caused by 
the QDC pedestal and broadened by dark current from the photomultipliers. 

To exclude events in the analysis of a detector channel where that channel was not hit by 
the electron beam, a cut threshold was dehned based on dedicated dark current measure¬ 
ments (i.e. runs without beam) taken over the course of the testbeam period. It was found 
that applying a threshold of 15 QDC bins above the QDC pedestal position eliminates at 
least 99.7% of the events in any dark current spectrum. For the analysis of the data, only 
events where selected in which both anodes of a channel yield a signal of at least 15 QDC 
bins above the mean of their respective pedestal. The part of the spectrum used in the 


analysis, compared to the spectrum before any cuts, is depicted in hgure 18b 


Figure 19 shows the signal of a single anode in both data and simulation. To convert 
the number of detected photons in the simulation to QDC bins, a gain of 1.3 ■ 10® (c.f. 


section 4.5) and the same digitiser resolution as in the data were assumed. 

In the data, the QDC signal has a noticeable tail to higher charges. The most likely 
explanation is the presence of secondary electrons in addition to the main beam electron, 
which would produce some additional light when they hit one of the quartz channels. 
It is currently not clear where these stem from. The wall thickness of the light-tight 
box surrounding the detector is only 4 mm, which does not account for enough shower 
electrons to explain the amount of additional light observed. Other possible contributions 
along the beam line have not been investigated so far. 
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QDC bins QDC bins (pedestal subtracted) 

(a) (b) 

Figure 18: Pedestal correction. Q An example for an uncorrected QDC spectrum is 
shown in light blue, recorded in anode 5 with the beam centred on channel 
A and a beam incidence angle of 45°. A dark current spectrum for the same 
anode is drawn in dark green (dashed line). (|^ The same QDC spectrum 
with shifted a:-axis to correct for the pedestal contribution. The light blue 
histogram contains all events, while the violet histogram (dashed line) is 
made up of only the events used in the analysis, i.e. events in which both 
anodes of a channel see a signal of at least 15 QDC bins above the pedestal. 


To account for the tail in the distribution, the QDC signals were htted with a £t function 
consisting of a Landau function convoluted with a Gaussian. The free parameters of the 
function are the width and most probable value (MPV) of the Landau function, the sigma 
of the convoluted Gaussian and a normalisation constant for the peak area. Whenever the 
testbeam data is compared with the simulation, the same £t is applied to the simulation. 
For the simulation, the width of the Gaussian is dominant, and the relation between MPV 
and the Gaussian width is consistent with the expectations from the statistical modelling 
used for the fast simulation introduced in section In the data, the Gaussian width 
is roughly twice as large as one would expect from noise in the amplihcation process 
alone. Since this broadening was not observed in a separate setup in the laboratory, it 
is attributed to the experimental environment at the testbeam, such as the long signal 
cables. The width of the Landau part of the £t function required to describe the large 
tails in the data is comparable to the Gaussian width. 


The most probable value of the htted Landau function was used as a measure for the 
amount of light in a channel, while the area parameter is related to the number of times a 
channel was hit by the beam. To illustrate the stability of the detector response and the 
fit procedure’s outcome, figure |20a shows the htted most probable value for both anodes 
of a channel for a number of data taking runs at identical beam position over a period of 
several hours. 


4.2 Horizontal alignment and beam profile 

To determine how well the prototype detector can be aligned from data alone, the detector 
position was scanned by varying the position of the movable table in x-direction in 1 mm 
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Figure 19: Fit of a Landau function convoluted with a Gaussian to the signal for Q 
anode 2 in a data taking with the beam centred on channel B, under an 
incidence angle of 30° and (|^ the corresponding simulation. The cut-off at 
low QDC bins is due to the pedestal subtraction cut. 
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Figure 20: (|aj) Fit results anodes 1 and 2 for 32 for consecutive measurement runs taken 
over a period of 11 hours (beam on channel B, under an angle of 30°). The 
observed light yield stays constant. B Landau area of anode 1 observed at 
different x-positions of the detector, htted with a box function convoluted 
with a Gaussian (as in equation 11). 


steps. Studying the observed number of events in each detector channel for the different 
x-positions allows to determine the position of the channel with respect to the beam, as 
well as derive information on the size and shape of the beam spot. For a pointlike beam 
and perfect alignment, the channel response would be expected to have a “box-like “ shape, 
whereas a hnite beam-size and tilts around the vertical axis will smear out the edges of 
this step response function. To determine the channel centres, the measurements for each 
anode were htted with a box function convoluted with a Gaussian. The box function is 
1 for all X within ±2.5 mm of its central value, to describe the prototype channel width 
of 5 mm, and the Gaussian is used to account for the smearing. The convolution is 
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channel 

anode 

channel centre 

anode 

channel centre 

A 

5 

(74.12 ±0.01) mm 

6 

(74.07 ± 0.01) mm 

B 

1 

(79.45 ±0.01) mm 

2 

(79.47 ±0.01) mm 

G 

4 

(85.54 ± 0.01) mm 

7 

(85.57 ±0.01) mm 

D 

0 

(91.10 ±0.03) mm 

3 

(91.15 ±0.02) mm 


Table 3: Channel positions fitted to x-scan data, with statistical errors from the fit. 

The central valnes for the two anodes per channel agree within 0.05 mm or 
better with each other. The positions are given in the coordinate system of the 
movable table’s control software. 


implemented as snm within fonr sigma of the Ganssian, i.e. 


P1+4-P2 

f{x)=Po- box(r,pi) 

r=pi-4-p2 


^ j 2-pi 


( 11 ) 


where the free parameters are a normalisation factor po, the centre position pi of the 
channel box fnnction, and the sigma p 2 of the Ganssian. As an example, the fit to one of 
the anodes is shown in fignre 20b The channel centre positions derived from these fits 


are listed in table The fitted centre positions to the data from the npper anodes of 
each channel all agree within 0.05 mm with the fitted positions for the respective lower 
anodes, indicating that the detector position can be determined within a precision of the 
same order, which meets the reqnirements for the horizontal alignment precision. 


In the fit of eqnationj^to the x-scan data of all eight anodes, the sigma of the convolnted 
Gaussians were fonnd to be a = 4.5 ±0.1 mm. By simulating x-scans with different beam 
spots and comparing them to the measurement, it could be verified that a beam profile 
with this extension would cause the smearing observed in the testbeam data. Therefore, 
all comparisons of the testbeam data to simulations were performed for simulations with 
a beam profile of a = 4.5 mm. 


4.3 Vertical alignment 

At the beginning of the testbeam campaign, the position at which the detector was centred 
with respect to the beam in vertical direction was determined by eye to be at a table 
position of ?/ = 76 mm (in the coordinate system of the steering software for the movable 
table). To check this estimate, data was taken at different ?/-positions over a range of 
40 mm, with the beam centred between channel G and D. Once the beam moves far 
enough down from the vertical centre that it enters channel D through the end face 
rather than the narrow side face, the path length for light production and thus also the 
amount of detected light will decrease rapidly. The same is true for channel G when the 
beam entrance point is moved upwards high enough. This effect was used to determine 
the centre in ^/-direction by comparing the data to a simulated y-scan. Two parameters 
need to be adjusted in this comparison: the first is a scaling constant to account for the 
difference in light yield in data and simulation; the second is an offset in the y-position. 
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Figure 21: Signal at different ^/-positions of the detector Q observed at the testbeam 
for both anodes of channel D (markers). The solid lines are the hts of a 
simnlated j/-scan to the signal of both anodes per channel, with a shift of 
the y-axis as a free parameter. e Simulated j/-scans with beam prohle 
c^beamspot = 0.1mm (left) and abeamspot = 4.5mm (right). 


i.e. how far the y-axis of the simulation (with the detector centre at position 0 mm) has 
to be shifted to match the data (with y-positions given in the coordinate system of the 
table software). 


A separate fit was done for both channel C and D. In these fits, both anodes of the 
respective channel were compared to the simulated |/-scan, allowing individual scaling 
factors for each anode in addition to the common offset in ^/-direction. Figure 21a shows 
the data for one channel together with the scaled and shifted simulated y-scans. The £t to 
the data of channel C returns a vertical centre ^/-position of 72.5 mm ± 0.1 mm, the £t to 
channel D a y-position 80.0 mm ±0.2 mm. The poor agreement is most likely due to large 
beam spot size of a = 4.5 mm. Therefore the distinct edges of the channel edges, which 
would be the most discriminating factor in the alignment, are washed out. This effect is 
illustrated in hgure 21b , where simulated |/-scans for both a = 0.1mm and a = 4.5mm 
are shown. As another possible explanation a misalignment in the x-z-plane has been 
considered: The ground plate on which the detector box was standing was only affixed 
to the movable table on one side, which could allow the whole setup to be tilted around 
the z-axis. This would induce correlations between the x- and ^/-coordinate, which would 
affect channel C and D differently and could lead to different vertical centre positions. 
However, in the simulation small misalignments in az (u^ < 5°) were found to have 
negligible impact on the light yield, and a larger tilt would have been visible by eye. 


The fit results do nevertheless conhrm that the first estimate by eye of a centre at //table = 
76 mm was accurate within a few millimetres. Under conditions where the vertical spread 
of the incident particles is smaller, a better alignment precision is expected]^ 


^More details on this can be found in 
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Figure 22: Angular dependence of light yield. Q Signal for both anodes of channel B 
for two consecutive angle scans, the hrst one from 30° to 60° (solid markers), 
the second one in the opposite direction (open markers). The second scan is 
shifted by 1.15°. 0 Shape of the angular response function for an angle scan 
at I/table = 70 mm, compared to the best-matching simulations at a a =4° and 
cr« =6°. 


4.4 Angle scans 


The prototype detector was equipped with a stepping motor to change the angle be¬ 
tween the incident electron and the detector channels. The angle was calculated from 
the stepping motor’s internal step counter. At the edges of the allowed movement range, 
a mechanical end switch was located, which stopped the motor movement when it was 
pressed. In some cases the stepping motor experienced difficulties on leaving the switch, 
causing a disagreement between the angle calculated from the step counter and the actual 
angle. This becomes evident by comparing the detector response for individual angle 
scans taken under otherwise unchanged conditions. Figure 22a shows the signal for both 
anodes of channel B for two consecutive angle scans, the hrst one from 30° to 60° (where 
the end switch was reached), immediately followed the second angle scan back in the 
opposite direction. In order to reach an agreement between the detector responses, the 
angle coordinate of the second scan had to be shifted by 1.15°. Three such bi-directional 
data sets were taken during the testbeam campaign. The angle correction required to 
obtain the best agreement between opening and closing angle scans for these three data 
sets were 1.15°, 0.69° and 0.78°, leading to the conclusion that the given angles for the 
testbeam data are afflicted with an uncertainty of 0(1°). 


The data of different angle scans taken during the testbeam campaign were compared to 
a number of simulated angle scans with different choices for the microfacet distribution 
parameter a a-, describing the surface smoothness of the quartz channels, and for the offset 
in the vertical direction Ay, to study the possibility to determine these parameters from 
the data: The surface smoothness of the quartz bars used in the prototype detector was 
not known. For a more polished surface, i.e. a smaller cTq,, more directed reflections occur. 
Consequently, a change in the incidence angle of the beam is expected to cause larger 
changes compared to reflections on a rougher surface. The vertical alignment within an 
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uncertainty of several millimetres also had to be taken into account. An offset in the 
vertical direction moves the central point of the beam crossing by and therefore 

^ o j sin a 

affects the detector response especially at small angles. 


In order to compare the shape of the detector response in the testbeam data to the 
different simulations independent of the absolute light yield, the response functions for 
both data and simulation were normalised to their respective signal at 45°, allowing for 
a small offset in angle to account for the observed uncertainty in the angle reading for 
the data discussed above. The current simulation does not describe the angular response 
satisfactorily for angles > 55°. To illustrate this, hgure 22b shows simulated angle scans 


with both cTq = 4° and aa = 6° compared to the data of an angle scan at |/tabie = 70 mm. 
The differences in the shape of the response function hint that additional factors must be 
present which are not well modelled in the current simulation. 

No dehnite conclusions could be drawn for the vertical offset, since no single value was 
found to describe the data of several angle scans signihcantly better than any other. For 
any given vertical offset, the best agreement was found for cTq, = 5° ± 1°. This was the 
case for all angle scans which were compared to the simulation, leading to the conclusion 
that a surface smoothness in this range would be a good description of the quartz bars 
used in the prototype detector. 


4.5 Light yield 


In order to compare the light yield observed at the testbeam to the predictions from 
the simulation, the gain of each photomultiplier anode has to be known. According to 
the “typical values “ from the datasheets, the gain for the R7600U-03 photomultipliers 
should be twice as high as for the R7600U-04 photomultipliers. In the data however, no 
discernible difference between the signal of these two photomultiplier types used at the 
testbeam was found. Therefore, rather than using the datasheet values, a separate test 
setup was used to determine the actual gain values of the PMTs used with the prototype: 
inside a light-tight box, a small amount of UV light was produced with an LED and further 
reduced by a hlter mounted in front of the photocathode. A more detailed description 
of the LED test setup can be found in 14 . The hlter strength and LED voltage were 


chosen such that ~95 % of the time only the dark current signal was observed.Assuming 
that the photons reaching the photomultiplier are Poissonian distributed, and that the 
quantum efficiency can be approximated by a binomial distribution with 20 % probability 
for detection, this means that in the ~90 % of the remaining 5 % of the events, in which a 
signal was observed, this signal stems from a single photon. After subtracting the pedestal, 
the remaining signal was htted with a Gaussian to determine the gain. The htted mean 
provides the mean number of electrons measured in the QDC, which corresponds to the 
gain in the single photon case. Data was taken at PMT supply voltages from 700 V to 
900 V. The determined gain followed a logarithmic behaviour as expected. As an example, 
the htted gain for all four anodes of one of the R7600U-04 photomultipliers is shown in 


hgure 23a 


Once the gain of the photomultipliers is known, the measured QDC signals can be com¬ 
pared to the simulation’s predictions for the number of detected photons. For each de- 
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Figure 23: (j^ shows the htted gain of one of the R7600U04 photomnltipliers for different 
supply voltages, with the dashed line indicating the typical values provided 
in the datasheet. B Ratio between the light yield calculated from the data 
(using the measured gain) and the number of detected photons predicted by 
the simulation, for all anodes of both photomultipliers, for measurements at 
three different incident angles. 


tector channel, data taken under three different angles (35°, 45°, 55°) with both pho¬ 
tomultiplier types was used for the comparison. For the conhguration with the highest 
light yield, at 35°, summing up the signal for both anodes of a channel yields 107 ± 12 
photoelectrons, while the simulation predicts ~70 photoelectrons for this angle. 

The ratio between the observed number of detected photons in the data and the predic¬ 
tions from the simulation, i.e. configurations are shown in 

’ #detected photonSgi„^^l,,yon ’ ^ 


figure 23b The alignment of the prototype detector at the testbeam with the precisions 


identified in the previous sections introduces systematic uncertainties in this ratio. While 
the error on the horizontal alignment is negligible, a systematic uncertainty on the light 
yield prediction from the simulation of — = 2 % for the angular alignment with Aax = 1° 


Ap _ 

P ~ 

considered. 


and ^ = 10 % due to the uncertainty of 0(5 mm) for the vertical alignment have to be 


On average, the amount of detected light calculated from the data is higher than expected 
based on the simulation by a factor of 1.7 ± 0.4(stat.) ± 0.2(syst.)- This might be partly an 
effect of additional light provided by secondary electrons. In addition to this, a number of 
assumptions where made in the simulation, such as the roughness of the quartz surface, 
and the thickness of the optical grease layer. In regard of this, an agreement better than 
factor 2 between the data and the light yield in the simulation is considered quite satis¬ 
factory, especially since the simulation underestimates the amount of light and therefore 
provides a conservative estimate of the measurable light. 


5 Conclusions 

Precise polarimetry is crucial for precision measurements at future lepton colliders, such 
as the International Linear Collider. At the ILC, two Compton polarimeters per beam 
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aim for an accuracy of 5V jV = 0.25 %, limited by the linearity and the alignment of the 
baseline gas Cherenkov detectors. 

A novel concept for detecting the Compton-scattered electrons using quartz as Cherenkov 
medium has been developed. Due to its much larger light yield per incident electron, 
it allows to resolve individual peaks in the measured Cherenkov light corresponding to 
different numbers of Compton electrons. This concept can improve the systematic uncer¬ 
tainty of the polarisation measurement substantially. In particular, it allows to control 
non-linearities in the response of the photodetectors used to detect the Cherenkov light, 
but also aids in alignment of the detector array. 

A four-channel prototype detector has been built and operated in a hrst testbeam cam¬ 
paign. The detector response to single electrons has been measured and compared to the 
simulation. The light yield per electron, which was a factor of 1.7 ± 0.4(stat.) ± 0.2(sygt.) 
higher than expected from the simulation, was found to be suitable for the construction 
of a detector capable of resolving individual peaks in the required dynamic range at the 
upstream polarimeter of the International Linear Collider. 

The testbeam data allowed a horizontal alignment of the prototype detector on the order 
of 0.05 mm. While the vertical alignment and study of the angular dependence of the 
detector would beneht from further investigation under better constrained conditions, 
i.e. with a smaller beam spot size and a more accurate determination of the vertical 
alignment independent of the data, they were sufficient to allow to determine that the 
surface roughness of the quartz used for the prototype is best described by a microfacet 
parameter (Tq, = 5° ± 1°. The noise level of the electronic setup at the testbeam was not 
adequate, leading to a signal width factor 2 larger than considered attainable. 

The testbeam results do not change the conclusions from the simulation studies signif¬ 
icantly. The differences between simulation and testbeam data in the angle scans are 
too small to affect the design optimisation. The larger observed light yield leads to an 
expectation of a possibly higher light yield for all detector geometries, which might allow 
to use smaller quartz channels. However, the observed light yield should hrst be verihed 
in future tests of the current prototype with multiple electrons and less noisy electronics 
to establish that the achieved light yield does indeed allow to resolve individual Compton 
electrons. 

Simulation studies extrapolating the testbeam results to ILC conditions showed that 
single-peak htting can compensate photodetector non-linearities up to 4% and allows 
a horizontal alignment to about 10 pm. The impact of tilts in the horizontal plane 
(0.07%/mrad) stays similar as for the baseline gas detector concept. The non-linearity 
compensation and improved alignment would reduce the effect of the two leading sources 
of systematic uncertainties on the polarimeter measurements signihcantly, to a level com¬ 
parable to the knowledge of the laser polarisation. This could reduce the total systematic 
uncertainty to even 5V jV = 0.20%, at least in case of the upstream polarimeter. In case 
of the downstream polarimeter, further studies are needed to investigate whether the low 
Cherenkov threshold of quartz poses a problem in the harsh background conditions. 
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